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Visualizing the elongated vortices in 7-Ga nanostrips 
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We study the magnetic response of superconducting y-Ga via low temperature scanning tunneling 
microscopy and spectroscopy. The magnetic vortex cores rely substantially on the Ga geometry, and 
exhibit an unexpectedly-large axial elongation with aspect ratio up to 40 in rectangular Ga nano¬ 
strips (width / < 100 nm). This is in stark contrast with the isotropic circular vortex core in a larger 
round-shaped Ga island. We suggest that the unusual elongated vortices in Ga nanostrips originate 
from geometric confinement effect probably via the strong repulsive interaction between the vortices 
and Meissner screening currents at the sample edge. Our finding provides novel conceptual insights 
into the geometrical confinement effect on magnetic vortices and forms the basis for the technological 
applications of superconductors. 

PACS numbers: 68.37.Ef, 74.25.Op, 74.78.Na, 74.25.Ha 


Magnetic vortices in mesoscopic superconductors, in¬ 
cluding symmetric disks, triangles, squares and rectan¬ 
gles m, have recently stimulated tremendous interest 
due to their fundamental importance in controlling su¬ 
perconducting energy dissipation, a significant feature for 
superconductor-based nanotechnologies. When a super¬ 
conductor is shrunk to nanoscale, exotic vortex shapes 
(i.e. giant vortex and antivortex-vortex molecule) and 
configurations may occur, but so far mostly inferred from 
indirect measurements [5] . Direct real space visualization 
of vortices in nanostructured superconductors remains 
challenging and rarely explored. A few scanning tunnel¬ 
ing microscopy (STM) studies present the preliminary 
evidence of giant vortex in Pb islands [6]-[8], but the ill- 
defined geometry of the samples investigated renders the 
issue quite complicated and difficult for a direct compar¬ 
ison with theory. 

Among all these superconducting nanostructures, rect¬ 
angular nanostrips or nanowires are especially appealing 
for superconducting electrical circuit and quantum inter¬ 
ference device. Regular oscillations of resistance and crit¬ 
ical current as a function of magnetic field have been re¬ 
spectively reported in superconducting InO and A1 strips 
laiini, which can be accounted for by the so-called Weber 
blockade theory m- Here each oscillation in either mag¬ 
netoresistance or critical current corresponds to adding 
a single vortex into the strip, in analogy with single elec¬ 
tron transport through quantum dots in the Coulomb 
blockade regime. Furthermore, the penetrating vortices 
are anticipated to configure into a one-dimensional chain 
at the center of strip noHis, because of the lower vor¬ 
tex potential energy there. Upon increasing field, the 
vortex chain may split into two or more buckled parallel 
rows HMH]. However, so far as we know, such vortex 
chains have little been directly visualized, and whether 


the model applies for the extremely narrow supercon¬ 
ducting nanostrips remains unjustified. 

Herein we report the direct visualization of unexpect¬ 
edly elongated vortex cores in superconducting y-Ga 
nanostrips by STM, and show that they are caused by 
the strong repulsive interaction between the vortices and 
Meissner screening currents at the sample edge. All ex¬ 
periments were conducted in a Unisoku ultrahigh vac¬ 
uum STM system equipped with molecular beam epi¬ 
taxy (MBE) for in-situ sample preparation. The MBE 
growth of superconducting Ga nanoislands has been de¬ 
scribed in detail elsewhere [T], and in the Supplemental 
Material HZ]. The measurements were carried out at 4.4 
K unless otherwise noted. A magnetic field up to 7 T 
can be applied perpendicular to the sample surface. The 
differential conductance dl/dV spectra and maps were 
measured by disrupting the feedback circuit at the set- 
point of U = 10 inV and I = 100 pA, sweeping the 
sample voltage Ug, and extracting the conductance using 
a standard lock-in technique with a bias modulation of 
0.3 mV at 987.5 Hz. 

Eigure 1(a) shows the constant-current topographic 
image of a rectangular Ga nanostrip, which has a length 
I ^ 440 nm and a width re ^ 64 nm. Note that the nanos¬ 
trip exhibits an almost perfect rectangle, which can easily 
be modeled in theory. The atomically resolved STM im¬ 
age on the nanostrip, shown in Eig. 1(b), reveals defect- 
free surface. The unit cell, marked by the blue rhombus, 
consists of two bright spots with a periodicity of uq = 
0.86 nm and an intersection angle of 76°, which are con¬ 
sistent with 7-Ga(001) plane [Eig. 1(c)] [18]. Eigure 1(d) 
plots a series of differential conductance dl/dV spectra 
acquired on the Ga nanostrip at various temperatures, 
normalized to the normal-state one above Tc (10 K). A 
superconducting gap with clear coherence peaks is visi- 
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FIG. 1. (color online) (a) STM topography of a rectangular 
Ga nanostrip with a height of 5.85 nm (45 ML) {Vs = 3.5 V, 
/ = 50 pA, 483 nm x 161 nm). Line profile taken along the 
red curve characterizes the height and width w of Ga nanos¬ 
trip. Both labels a and b denote respectively the orientations 
of the short and long side of the rectangular strips through¬ 
out this paper. Dotted rectangle marks the region where all 
ZBG maps are taken in Fig. 2(a). (b) Atomic resolution im¬ 
age (l/g = 20 mV, / = 50 pA, 10 nm x 10 nm) of the Ga 
nanostrip. The bright spots correspond to the Ga atoms in 
the top layer, (c) Schematic crystal structure of y-Ga. Along 
the c(OOl) axis, each unit cell consists of four atomic layers. 
For clarity, the spheres in each layer are coded in different 
opacity. The brightest Ga spheres at the top layer form a 
dimmer-like structure, consistent with the STM image in (b). 

(d) Temperature-dependent dl/dV spectra on Ga nanostrip 
(black dots) and their best fits to the Dynes’s expression (red 
curves), with the broadening parameter F of 0.1 ~ 0.2 meV. 

(e) The extracted gap magnitude A versus temperature (black 
squares). The red line shows the fit to BGS gap function. 


ble at 2.5 K. At elevated temperature, the gap gradually 
reduces aud both the cohereuce peaks are suppressed, as 
expected. The supercouductiug gaps cau be fitted by us- 
iug Dyues’s expressiou m, as illustrated iu Fig. 1(d). 
Figure 1(e) plots the extracted supercouductiug euergy 
gap A at various temperatures. Fittiug the data to BCS 
gap fuuctiou yields Tc = 6.80 K aud A(0) = 1.28 meV. 

Haviug ideutified the lattice structure aud supercou- 
ductivity, we theu focus ou the maguetic respouse of 
7 -Ga uauostrip. To visualize the maguetic vortices iu 
real space, we map the spatial variatiou of zero bias cou- 
ductauce (ZBC) ou the dotted rectaugle iu Fig. 1(a), a 
commou approach for STM imagiug of vortices [OHHl HDJ- 
[27]. This techuique takes advautage of the ZBC coutrast 
withiu aud outside the maguetic vortex cores, aud has a 
higher spatial resolutiou (^ ^). Here, the suppressed su- 
percouductivity withiu the vortex cores leads to a higher 
ZBC value, leaviug supercouductivity ouly outside the 
vortex cores with lower ZBC. Thus the regious with eu- 
hauced ZBC iudicate the emergeuce of maguetic vortices. 
Figure 2(a) illustrates such maps at various fields raugiug 
from 0.3 T to 1 T. Two iutriguiug pheuomeua are imme¬ 


diately uoticed. The first pheuomeuou, which is also the 
more iuterestiug aud uuexpected oue, is that the mag¬ 
uetic vortex cores expaud aud compress aloug the louger 
{a axis) aud shorter {b axis) sides of the rectaugular Ga 
uauostrip, respectively, regardless of the field. This cou- 
sequeutly leads to the auisotropic vortex cores with great 
axial elougatiou. The elougatiou ratio r] is estimated to 
be 15 at 0.435 T, aud iucreases further with the maguetic 
field. At 1.0 T, mauy vortices merge together aud the su¬ 
percouductivity is almost completely killed, pushiug the 
whole islaud iuto the uormal state. Figure 2(b) reveals 
the field-depeudeut ZBC profiles across the vortex cores. 
Aloug the b axis, the existeuce of ouly oue extremum iu 
every ZBC profile excludes reasouably the possibility of 
vortex chaiu for a siugle iudividual yellow stripe. Ou the 
other haud, ZBC oscillates aloug the a axis, with each 
peaked regiou represeutiug siugle or multiple (e.g. giaut 
vortex) of the maguetic flux quautum ^o- The uumber 
of peaks iuvariably appears smaller thau, but teuds to 
close the applied flux the field iucreases. Above 

0.7 T (^/^o = 10T), the ratio betweeu them becomes 
greater thau 0.5, ruliug out the possibility of giaut vortex. 
We therefore suggest that each peak (or yellow stripe) 
actually deuotes a siugle flux quautum uamely the 
vorticity L = 1 [Fig. SI]. 

Secoudly, we fiud that the maguetic vortex first peu- 
etrates iuto the Ga uauostrip iu a rather high filed, i.e. 
0.435 T iu Fig. 2(a). This is caused by the fiuite eu¬ 
ergy barrier for vortex eutry, cousistiug of the geometric 
[28] aud Beau-Liviugstou (BL) surface barrier [29]. Both 
barriers couspire to preclude the peuetratiou of vortices 
at low field. The first vortex peuetratiou occurs ouly 
wheu the surface supercouductiug curreuts exceed the 
pair-breakiug curreut, locally weakeuiug supercouductiv¬ 
ity aud theu allowiug the uucleatiou of a vortex. 

The auisotropic iuterual vortex structure coustitutes 
the major fiudiug iu this study: their large elouga¬ 
tiou aud the mechauism behiud are uuprecedeuted. For 
auisotropic vortex cores, such as sixfold star-shape [20] . 
fourfold [23H25] aud twofold symmetry [211 EH EO] , to 
uame a few, seeu previously iu certaiu uucouveutioual 
supercouductors, the auisotropy iu either the Fermi sur¬ 
face or the supercouductiug gap has beeu suggested to be 
the primary cause. However, this is uot the case here: iu 
7 -Ga the supercouductiug gap cau be well fitted by the 
isotropic BCS theory [Fig. 1(e)], aud uo gap auisotropy 
is iuvolved. To shed more iusight iuto the possible mech¬ 
auism behiud, we have iuvestigated the vortices iu more 
Ga uauostructures, aud fiud that the vortices iuvariably 
elougate ouce if the width of Ga uauostrips is reduced to 
below ^ 100 um. Iu coutrast, usual isotropic circular vor¬ 
tices develop iu a roughly rouud-shaped Ga islaud with a 
larger lateral dimeusiou [Fig. 3(a)], which although shows 
the same lattice structure aud supercouductivity as Ga 
uauostrips. Showu iu Figs. 3(b-i) are the isotropic vortex 
cores, which bear stroug similarities with those observed 
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FIG. 2. (color online) (a) Normalized ZBC maps at various magnetic fields, acquired in a field of view of 380 nm x 51 nm on 
the Ga nanostrip of Fig. 1(a). The normalization was performed by dividing these maps by the normal-state ZBG value at the 
vortex core. Yellow regions with enhanced ZBG correspond to the vortex cores, all showing enormous elongation along the b 
axis. The small ZBG enhancement on the right side Ga nanostrip at 0.435 T probably indicates that we are probing a critical 
state justly before the penetration for the second vortex, (b) ZBG profiles across vortex centers in (a). The ZBG presents 
oscillations along the a axis, while it shows a single extremum along the b axis. Here 4>/4>o and L denote the applied magnetic 
flux and total vorticity in the strip, respectively. 



FIG. 3. (color online) (a) STM topography (14 = 3.0 V, / = 50 pA) of a larger Ga island (350 nm x 280 nm) with a height 
of 3.39 nm ( 26 ML). The dotted square marks the region over which all ZBG maps in (b-i) are acquired, (b-i) Magnetic Field 
dependence of ZBG maps, acquired in a field of view of 140 nm x 140 nm in (a). The isotropic vortex cores are unchanged 
respective of the fields, (b) ZBG map (140 nm x 140 nm) taken on a round shaped Ga island with a height of 3.39 nm (26 
ML), showing a highly isotropic vortex core, (j) Radial dependence of a{r) across the isotropic vortex, which are color-coded 
to match ZBG map in (b). Here a(r) has been normalized to optimize their fit to the inserted expression by minimizing the 
residual, as shown by the red curve, (k) Magnetic vortex-induced ZBG enhancement on various Ga nanostructures. 


in conventional superconductors such as Pb [22]. This 
completely excludes the anisotropy in intrinsic electronic 
structure of y-Ga as a possible cause for the elongated 
vortices. Instead, the vortex elongates mainly via a geo¬ 
metrical confinement effect. 

Figure 3(j) plots the radial dependence of vortex- 
induced ZBC [or cr(r)]. Based on the Ginzburg-Landau 


expression for the superconducting order parameter A 
near the interface between a superconductor and a nor¬ 
mal metal, the radial ZBC profile across the vortex core 
should obey [21] 

<^{r) = cToo + (1 - croo)(l - , (1) 
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FIG. 4. (color online) (a) Schematic of vortex axial elonga¬ 
tion. Je and Js denote the Meissner screening currents at 
the sample edge and around the vortex, respectively. Black 
thick arrows label the repulsive force between the vortex and 
Meissner screening currents Je. (b) Magnetic held and 7 -Ga 
geometry dependence of vortex elongation 77 . The extracted 77 
for l/w — 1.0, 2.1, 6.9 and 8.5 are coded magenta, red, black 
and blue, respectively. 


core cannot be ignored. Moreover, the nanostrip thick¬ 
ness (several nanometers) is considerably smaller than 
the penetration depth A, vortices are of the Pearl rather 
than Abrikosov type [32]. The vortex interactions with 
the Meissner screening currents at the sample edge (Je) 
becomes long ranged and prominent. In terms of the 
counterpropagation between Je and the screening current 
Js encircling the vortex (left panel in Fig. 4(a)), the inter¬ 
actions are dominantly repulsive and pronounced along 
the a axis. To minimize this repulsive interaction, the 
vortex current Jg will deform or more specihcally elon¬ 
gate along the b axis until they are fully exploring the 
sample geometry (right panel in Fig. 4(a), less energeti¬ 
cally costly), where the Je-Js repulsive interaction along 
the b axis gets stronger as well. Consequently, an equi¬ 
librium with the strong deformation in Jg and thus ZBC 
or vortex cores is reached, as observed. 


where cr(r) is the normalized ZBC away from the vortex 
core and r is the radial distance from the vortex center. 
A best fit of ZBC profile cr(r) to Eq. ( 1 ) yields a super¬ 
conducting coherence length ^ = 24.6 ± 1.1 nm at 4.4 
K. Figure 3(k) further qualifies the vortex-induced ZBC 
enhancement as a function of the vorticity L in the round 
Ga island and rectangular Ga nanostrips, seen more vi¬ 
sually in Fig. S 2 . The linear relationship between them, 
guided by the dashed line, indicates the equivalence (ex¬ 
cept for the shape) among all observed magnetic vortices, 
regardless of the geometry. 

We now consider possible explanations for the elon¬ 
gated vortices in the Ga nanostrips. Firstly, as the width 
w of Ga nanostrips is comparable to their Fermi wave 
length Ap, quantum confinement will become important 
and might lead to sizable spatial modulation in A across 
the width of nanostrips |3|, forming a single or a multi¬ 
ple of weak superconducting thin slices normal to the a 
axis. As a result, vortices, as locally quenched supercon¬ 
ductivity, will most probably reside inside the weak su¬ 
perconducting slices (more energetically favorable), and 
may enormously elongate along the slices EZIEI]. How¬ 
ever, our Bogoliubov de-Gennues (BdG) self-consistent 
calculations show that this explanation from the quan¬ 
tum confinement effect is applicable only if the Fermi 
wave length Ap is larger than about eight times of the 
lattice constant. In a usual metal, however Ap is com¬ 
monly comparable to its lattice constant (e.g. Ap r\j 0.4ao 
in 7 -Ga) [1], as indicated by considering free electron ap¬ 
proximation. Therefore, the quantum confinement effect 
might not play the primary role in the vortex elongation 
observed here. 

Alternatively, another possible origin of elongated vor¬ 
tices is from the geometric confinement effect imposed 
by narrowness of the rectangular Ga nanostrips. In the 
Ga nanostrips studied here whose widths are of the or¬ 
der of ^ (e.g. w ^ 2.5^ in Fig. 2), the finite size of vortex 


From considering Je-Js repulsive interactions, one may 
further argue that vortex elongation 77 could enhance as 
the ratio l/w and magnetic field increases, because Je in¬ 
creases with increasing field. In addition, once if multiple 
vortices enter into the nanostrips, the strong interaction 
between vortices will not only compress a certain vortex 
along the a axis as well and contributes to enhance 77 , 
but also configure all vortices parallel to the b axis, as 
observed. Figure 4(b) summarizes the vortex elongation 
77 as function of the field B and l/w oi the 7 -Ga nanos- 
tructrues. Clearly, 77 increases abruptly with B and l/w. 
For example, in the Ga nanostrip with the aspect ratio 
l/w of 8.5, 77 is as high as 40.5 at 0.7 T. Therefore, even 
though further detailed theoretical investigations will be 
needed to understand the vortex elongation in a quanti¬ 
tative way, our experimental findings support that the re¬ 
pulsive interaction between the vortices and sample edge 
screening current is a reasonable explanation of the vor¬ 
tex elongation in Ga nanostrips. 

To summarize, our real space visualization of unusu¬ 
ally elongated vortices in 7 -Ga nanostrips sheds light on 
Pearl vortex dynamics in the extremely thin and narrow 
superconducting nanostrips {w < 4^), where the geomet¬ 
rical confinement effect may deform the internal vortex 
core structure and then configure vortices in a brand- 
new way. The elongated vortices observed here provide 
a novel platform to study exotic vortices in nanoscale su¬ 
perconductor. From the viewpoint of superconductor ap¬ 
plications, the strong geometrical confinement could im¬ 
pose deep potential well for the vortex transverse motion 
and immobilize the vortices, providing a promising strat¬ 
egy to design nanoscale superconducting devices working 
at high magnetic field, a long-standing dream of super¬ 
conducting research. 

This work was supported by National Science Founda¬ 
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SUPPLEMENTAL MATERIAL 

Commercially purchased Si(lll) wafers with a resis¬ 
tivity of 0.01 0-cm were used as substrate, and the clean 
Si(lll)-7 X 7 surfaces were obtained by flash annealing 
to ^ 1200°C while keeping the vacuum better than 1 x 
10-9 Xorr. We then prepared Ga nanoislands by evap¬ 
orating Ga (99.999%) sources from a standard Knudsen 
cell at a nominal beam flux of approximately 0.4 mono- 
layer (ML, ^ 0.13 nm)/min, during which the Si(lll)- 
7x7 substrate was kept at 110 K. Subsequent post¬ 
annealing at a slightly higher temperature of 150 K led 
to many Ga nanostructures with varying geometries and 
the lateral dimensions [T]. Analogous to MBE growth of 
Pb on Si(lll)-7 x 7 [5], the substrate just beneath Ga 
islands consists of a slightly disordered Ga wetting layer 
(1^2 ML), which exhibits no superconductivity down to 
the minimum measurement temperature of 2.2 K. 
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field B for three Ga nanaoislands with varying geometry and 
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Figure SI depicts the evolution of vorticity L (the num¬ 
ber of individual isolated stripy or round yellow features 
in ZBC maps) with varying magnetic field B. Above a 
threshold field 5^, depending on A or the width W of 
superconducting strips, L increases almost linearly with 
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as anticipated for the gradual penetration of mag¬ 
netic vortices into Ga islands [3] . By extracting the slope 
(dashed lines) and island area A, we estimate that each of 
such features accommodate the magnetic flux of 1.18^o 
(left triangle), 1.25^o (upper triangle) and 1.47<ho (cir¬ 
cle) , respectively, close to a single magnetic flux quantum 
4>o and never exceeding 24>o. This strongly supports a 
single magnetic vortex nature of the observed stripy or 
round features in Figs. 2 and 3. The small discrepancy 
from 4>o most possibly originates from an overestimate 
of the effective island area A for nanoscale superconduc¬ 
tors. Compared to Nb strips previously reported [3], 
appears larger in Ga strips we studied here, due to their 
extremely small dimensions. 
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FIG. S2. Differential conductance dl/dV spectra straddling 
(a) a round magnetic vortex core (Fig. 3, 0.15 T) and (b) an 
elongated vortex core (Fig. 2, 0.435 T), respectively. Inserted 
are the two vortices, with the dashed lines indicting where 
the dl/dV spectra are acquired. As anticipated, the super¬ 
conductivity is significantly suppressed with enhanced ZBC 
when approaching the vortex cores, regardless of their inter¬ 
nal core structure.The small F^-near conductance reduction 
in (b) is primarily caused by a large (~ 12 nm) spatial in¬ 
terval between the neighboring dl/dV spectra, which fails to 
acquire the spectrum justly at the vortex core center. 














































